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ABSTRACT 

The long exposure times of the HST Ultra-Deep Field plus the use of an empirically 
derived position-dependent PSF, have enabled us to measure a cardioid/displacement distor- 
tion map coefficient as well as improving upon the sextupole map coefficient. We confirmed 
that curved background galaxies are clumped on the same angular scale as found in the HST 
Deep Field North. The new cardioid/displacement map coefficient is strongly correlated to 
a product of the sextupole and quadrupole coefficients. One would expect to see such a cor- 
relation from fits to background galaxies with quadrupole and sextupole moments. Events 
that depart from this correlation are expected to arise from map coefficient changes due to 
! lensing, and several galaxy subsets selected using this criteria are indeed clumped. 
Of 
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^ ■ 1. Introduction 



This letter expands upon previous work of two of us (Irwin and Shmakova 2003b, 2005), hereafter 
[IS] , aimed at identifying lensing signals created by small-scale structure and substructure, the ultimate 
purpose being the development of a unique small-scale cosmological signal. 

In [IS] we described and developed a non-linear lensing method based on a complex power series 
lensing distortion map. The map is derived from a symplectic coordinate transformation from a reference 
frame at the telescope to a reference frame at the source galaxy. It has a succinct complex power series 
form: w s = w + aw + bw 2 + d\ww + d^w 1 , where w = x + iy represents the position relative to the 
beam centroid. This map is a generalization of the usual linear weak lensing relationship. Each term 
in the series is a rotational eigenvector and creates a distinct distortion of an azimuthally symmetric 
source galaxy. There are two well-known linear terms: an amplification ("convergence") and a term 
that produces a quadrupolar distortion ("shear"). For lensing a(= a L ) = —7. 

The three 2nd-order map terms are: a sextupolar triangle-like distortion, a cardioid distortion, 
and a quadratically increasing displacement of circles. Since lensing simultaneously produces both the 
cardioid (g^) and displacement (di) terms, they can be denoted by a single coefficient (d), with d\ = 2d 
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and d>2 = d . We call the so-combined (i-terms the "cardioid/displacement" term. See fig. 1 for examples 
of these shape distortions. 1 
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Fig. 1. — Shape deformations coming from different terms of a map. The dashed circles correspond to 
an original isophote. These pictures demonstrate the deformation due to the presence of single terms: 
either a x = 0.3, b x = -0.06, d x = 0.02. 

The lensing distortion map embodies a distinction between intrinsic galaxy shapes and allowed 
lensing distortions of that shape. The power of this approach will become particularly evident here as 
we expand our consideration to include all 2nd-order map coefficients. The interrelationship of these 
coefficients allows one to infer whether a particular background galaxy is likely to have been lensed. 

In [IS] the lensing distortion map formalism was applied to the Hubble Deep Field North (HDFN). 
PSFs generated by ? were convolved with mapped radial profiles. Parameters for the map and the radial 
profile were simultaneously determined by a least-squares fit to the background galaxy intensity matrices. 
Curved background galaxies, defined by the relative orientation (for that galaxy) of the quadrupole and 
sextupole map coefficients, were found to be spatially clumped when compared to randomly selected 
background galaxy subsets, suggesting that some fraction of curved galaxies are the result of lensing by 
foreground haloes and/or subhaloes. We find it improbable that the observed clumping could arise from 
intrinsic alignment of background galaxies since the members of identified clumps have a wide range of 
z- values. 

We report here the results of extending that analysis to the HST Ultra-Deep Field (UDF). Due to 
the longer UDF exposure and access to an empirically derived time-averaged field-dependent PSF for 
the Hubble wide- field camera (WFC) (Anderson and King 2006), we have again observed the clump- 
ing of curved galaxies and we have successfully measured an additional second-order map coefficient 
representing the cardioid/displacement distortion. 

We find a strong correlation between the the cardioid/displacement coefficient and a product of 
the quadrupole and sextupole map coefficients. Such a correlation is expected from a fit to background 
galaxy shapes. Galaxies that do not comply with this correlation may be considered as candidates for 
quadrupole or sextupole lensing or both. We present results supporting this hypothesis. 



1 Other authors (Bacon et al. (2005); Goldberg, and Bacon (2005)) have referred to a related, though differently 
parameterized, concept as Flexion I. 
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2. UDF data processing 

The F606W UDF data set which we analyze in this paper consists of 112 half-orbit exposures. Each 
exposure was taken at a different pointing, with two primary orientations that were rotated by ~ 90°. 
We transformed each frame of the field into a distortion-free frame which one of us (Anderson 2005) 
had constructed for the WFC. 

Using the distortion-free transform of the first UDF frame as a frame of reference and bright 
point-like galaxies as locaters, we next determined linear transformations between the distortion-free 
transform of each of the remaining frames and this reference frame. These transformations were then 
used to map the center of each pixel in each frame into the reference frame. We removed the sky from 
the images by examining pixels in 800x800 patches and setting the background to zero. 

The construction of the super-image is an iterative process. The super-image is chosen to be 
2x super-sampled with respect to the image pixels. To construct this 8200x8200 super-image, we 
proceed pixel-by-pixel finding the closest corresponding pixel in each individual frame. This provides 
112 independent estimates of the flux for each pixel. Since the individual frames have slightly different 
exposure times, we divide the flux by the ratio of the image's exposure to a reference-image exposure 
time of 1200 seconds. We employ simple sigma-clipped averaging to find a representative value for that 
pixel, obtaining a first super-image estimate. Only one out of four pixels is statistically independent. 

So that the super image can be interpolated without high frequency noise we smooth the first 
composite image with the kernel 

"0.0625 0.1250 0.0625" 
0.1250 0.2500 0.1250 . 
.0.0625 0.1250 0.0625. 

Then for each pixel in each individual exposure, we find the difference between it and the interpolated 
super-image (from the first iteration), thereby obtaining 112 frames of residuals. We add the sigma- 
clipped average of these residuals to the original pixel value, to get the final value of each super-image 
pixel. 

To ensure that the galaxy shapes we are measuring are intrinsic to the galaxies and not artifacts of 
the detector, we propagate the model PSFs through the same stacking procedure as used for the galaxy 
images. To get a PSF appropriate for our super images, we generated an array of artificial stars in 
our reference frame. The stars had total fluxes of 10 5 photons and were separated by ~ 100 pixels and 
placed with pixel phases of (0,0), (0.5,0), (0,0.5) and (0.5,0.5). We used the inverse transformations 
to determine where each of these stars landed in each of the 112 exposures, and then used the spatially 
variable PSF model of Anderson and King (2006) to generate an artificial set of exposures with this 
grid of stars in them. 

Finally, we combine together the PSFs at different pixel phases to get a x4 representation of the 
PSF. It is this x4 PSF which we convolve with a x4 mapped radial profile to find the best fit to each 
galaxy image. 

Since the fit procedure to find map parameters is designed to analyze galaxies with a single dominant 
peak, we have developed software to identify and select such galaxies. This software identifies all peaks 
within each galaxy (previously selected using SExtractor softwareBertin and Arnouts (1996)) and based 
on criteria such as separation, flux, and pixel height, decides whether these peaks should be considered 
part of a single peak. After all peaks are processed then, based upon a second set of criteria such as total 
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flux and maximum-to-minimum height ratio, satisfactory central-peak regions are chosen to represent 
selected galaxies. 

The fit-method is the same as that described in Section 3.3 of [IS]. The radial profile is taken to be 
a quadratic polynomial in r 2 times a Gaussian. A condition is imposed that the polynomial be positive. 
The map has four pairs of parameters: the centroid, the real and imaginary quadrupole coefficients, the 
real and imaginary sextupole coefficients, and the real and imaginary cardioid/displacement coefficients. 
The procedure begins by determination of the full-width-half-maximum of the image and its ellipticity. 
Then, without yet introducing the PSF, the radial profile and map coefficients are fit using a least- 
square difference with the selected region of the galaxy intensity matrix. Convolutions with ever larger 
footprint PSFs are similarly fit until the final PSF is convolved with the mapped radial profile and fit. 

The distributions of the magnitudes for the quadrupole, sextupole, and cardioid/displacement 
coefficients of background galaxies (z > 0.5) are shown in fig. 2. The radial spreading caused by the 
PSF causes the map coefficients to be smaller then they are for the pre-PSF shape. 

p p p 
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Fig. 2. — Distributions of the magnitude of the quadrupole coefficients (a, left panel), the sextupole 
coefficients (b, center panel), and the cardioid/displacement coefficients (d, right panel), of background 
galaxies using the radial- fit method without PSF (dashed) or the model method with the PSF (solid). 



3. Clumping of curved galaxies 

We looked for evidence of a spatial clumping of both curved and aligned galaxies as found in our 
HDFN analysis. 2 The UDF was observed with the WFC camera which has a 2x smaller pixel size than 
the WFPC used for the HDFN. Because the total exposure time is about lOx longer, we were able to 
resolve 3 times more background galaxies per sq. min. We have used a different image stacking process 
and a different, empircally based PSF determination method. Strikingly, we find clumping of curved 
galaxies in the UDF at the same spatial angular scale as found in the HDFN. If lensing is the source 
of the clumping, it follows that the scattering lenses are similar in size and mass distribution. The 
probability of the observed clumping being random in the UDF is found to be ~ 1%, as compared to 
our result of 4% in the HDF. The distribution of the number of neighbors in circles of radius R=500 
pixels is shown in fig. 3. 1% is the probability that randomly chosen subsets (of all background galaxies) 



2 A "curved" ("aligned") galaxy was denned as having a minimum (maximum) of its sextupolar distortion within 10° 
of one of the two minima (maxima) of its quadrupolar distortion. 
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with the same number of galaxies as the curved set have as many galaxies with N=8 or more neighbors 
as the curved set. 

The clumping we observed for aligned galaxies in the HDFN is very weak in the UDF. 

p R = 500 pixels 
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Fig. 3. — The foreground bars show the distribution of the number of curved neighbors in circles of radius 
R=500 pixels for the set of all curved galaxies. The background shows the same distribution averaged 
over randomly chosen subsets of galaxies with the same number of members. 500 pixels corresponds to 
12 arcsec. Note the consistent shift to higher numbers of neighbors. 



4. Correlations in background galaxy shapes 

For our UDF images we have been able to measure the magnitude and orientation of an additional 
2nd-order lensing-distortion-map coefficient, which we have referred to as the cardioid/displacement 
coefficient, or d-term. The d-term distorts the variable r| = wsws in the argument of an azimuthally- 
symmetric light stream profile with an amplitude proportional to r 3 cos(9 — 9d). To our knowledge this is 
the first measurement of this coefficient with sufficient precision to exhibit the fundamental correlation 
we now describe. 

We find the orientation of the d-term coefficient (9d) to be correlated with the ah orientation (9b— 9 a ). 
Scatter plots of 9 d vs. (9b — 9 a ) are shown in the left and center panels of fig. 4. The center plot, in 
which the correlation is obvious, includes the PSF during the fit whereas the left plot, which barely has 
an indication of the correlation, ignores the PSF, illustrating the significance of the accuracy of the PSF 
treatment. We have studied the properties of the PSFs and find no evidence there for this correlation. 

A projected distribution of galaxy count vs. 9^ — (9b — 9 a ) is shown for all galaxies in the left 
panel of fig. 5, and for curved galaxies in the center panel. For comparison, in the right panel we show 
the shape that would be predicted for an underlying delta-function correlation that incorporates the 
noise as predicted from a Fisher matrix study of the determination of these angles. The widths are 
similar, suggesting that the underlying correlation, though not necessarily a delta function, is indeed 
considerably sharper. 

As the reference frame of the galaxy is rotated, the angles 9a and 9b — 9 a increase linearly, explaining 
why the correlation in the center panel of fig.4 lies along a line of unit slope, and making it evident 
that this correlation is invariant under rotation. A feature of the PSF or an effect such as CTE (charge 
transfer efficiency), would not possess rotational invariance, and thus could not lead to such a correlation. 

Another perspective on this correlation can be seen by plotting the real (imaginary) part of d versus 
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the real (imaginary) part of ab. This is done in the right panel of fig. 4. One sees a general tilt to this 
diagram along the line d = §a&. The factor 3/5 has been determined, to a few percent, by a least-squares 
fit of a line to the data. 

As a guide to find further correlations, one can use the moment methods of section 3.2 in [IS] to 
find the relationship d — ^ad = ~ab — |(1 — |a| 2 )m 2 i where m 2 i = and M nm = J w n w m idA. This 
equation was derived assuming an azimuthally symmetric root galaxy was mapped to fit a background 
shape. M21 and M22 are moments of the image of the root galaxy after the mapping. 

An equation of this form is to be expected from general considerations. The quantities d, b, and 77121 
are small enough that they occur only in first order. Furthermore as the reference coordinate system 
is rotated, d, db, m 2 i, a 2 b, ad, and am 2 i all rotate as vectors, so only these quantities can be involved. 
The terms with 777,21 and 077121 can be added together and treated as the unknown. And only one of the 
two quantities ad and a 2 b need be retained, as they are redundant. 

Minimizing the (combined) 777,21 term results in the approximate equation 

4-3 1 

d — -ad ps -db — -(1 — |a| 2 )7772i. (1) 
5 5 5 

The 3/5 ab and 4/5 ad are good fits to the data, and an additional a 2 b term, if present, has a very small 
coefficient. The 1 — \a\ 2 dependence of the coefficient of 77721 was adopted from the moment equation. 
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Fig. 4. — The angle of the d-term vs. the angle of a6 (#& — O ), without including the PSF (left 
panel) and including the PSF (center panel). The right plot shows Re[d] vs. Re[db]. 



5. Evidence for cardioid/displacement lensing 

The utility of the lensing-distortion-map parameterization is that if a background galaxy is subse- 
quently lensed, this correlation, expressed in terms of the observed total a ps a^ + a s and b ~ b^ + 65, 
where L indicates lensing and S denotes the fit to the background source galaxy, would break down 
even if the lensing is quadrupole lensing alone. In other words, we know that in the presence of lensing, 
the map coefficients representing the shape of the source galaxy will add linearly to the lensing map 
coefficients. Since a contribution from lensing can be expected to change these coefficients, lensing will 
disturb this correlation condition. This will show up as a larger magnitude of 77721 when it is calculated 
using eq.l, because the change in 3/5 ab will be vectorially added to the existing 777 2 i. 

Hence if we choose those galaxies that are in the tail of the |777 2 i| distribution (see the left panel 
of fig.6) they should be more likely to be lensed galaxies. We chose the cut at |777 2 i| > 0.023 to end up 
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Fig. 5. — The left (center) panel shows the distribution of the angular difference \9d — (9b — 9 a )\ mod 
180° for all (curved) galaxies. The right panel shows our estimate, based on Fisher matrix calculations 
for each galaxy, of the expected measurement spread in |6^ — (Oj, — 9 a )\ for a ^-function correlation. 

with 200 galaxies. Sets with less than 200 members have poorer statistics for clumping studies. The 
neighbor distribution for galaxies with |m 2 i| > 0.023 is shown in the central panel of fig. 6. It has a very 
small probability of occurring randomly, less than 1%. 

Additionally we have investigated the small but tantalizingly distinct hump in the distribution of 
the curved galaxies shown in fig. 5 near 6d — (6b — 9a) = 180°. For these galaxies the sign of the 
coefficient of d is opposite to that of ab. This is the orientation for lensing events that are non-zero for 
all three lensing coefficients, and for which the d-term coefficient would be responding to the presence 
of £' (the gradient of the projected mass density function) in the lensing plane. 

Since this is a small set of galaxies, it is difficult to convincingly establish lensing as a probable cause 
of this hump by a study of clumping probabilities. However we have done that, and when compared to 
randomly chosen curved galaxies of the same size, it has a probability of arising randomly of ~ 3%. It 
is perhaps more convincing to look directly at the neighbors distribution, as shown in the right panel 
of fig.6. The neighbors distribution is compared for three distinct subsets of curved galaxies having 
the same number of members. The galaxy subset with 150° < 9^ — (9b — 9 a ) < 180° clearly has larger 
numbers of neighbors. Though this is suggestive, we would not claim to have conclusively observed E' 
lensing. 



6. Conclusions 

We have confirmed in the UDF the spatial clumping of curved galaxies observed by [IS] in the HST 
Deep Field North. The clumping radius is the same as found in the HDFN with a smaller (now ~ 1%) 
probability of occuring randomly for the UDF. We are now confident that we have observed nonlinear 
lensing due to small-scale structure. 

We have measured, for the first time, the cardioid/displacement map coefficient (d-term) and its 
correlation with the quadrupole and sextupole. This correlation can be explained, in the absence of 
lensing, as arising from an M^ moment in the source galaxies, which gives rise to a fe-term and also, 
through ab, to a <i-term in the map. Galaxies not exhibiting this correlation are expected to be lensed, 
and we have presented evidence for the spatial clumping of the set of 200 galaxies having the largest 
values for the magnitude of 777.21. The probability of this clumping being random was < 1%. 
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Fig. 6. — The left panel shows the distribution of the magnitudes of 777.21 as calculated from eq.l. 

Large 1 777,21 1 is a signal for a violation of the <i-term correlation. The central panel shows the neighbor 
distribution for 225 galaxies with 1 777,21! > 0.023. The background shows the neighbor distribution 
averaged over randomly chosen subsets of galaxies with 225 members. Note the consistent shift to 
higher numbers of neighbors. The probability of this shift occurring randomly is < 1%. The right 
panel shows the distribution of the number of neighbors for three subsets of the curved galaxy 

distribution with 90 members: 1) — 4.3° (background), 2)45 — 100°, and 3)100 — 180° (foreground). 

We have also used the orientation of the d-tevm with respect to the direction of ab to identify a set 
of about 80 galaxies that could have been lensed by the density derivative, As would be expected 
from strongly lensed galaxies, this subset consists almost exclusively of curved galaxies and appears to 
be spatially clumped. 

As larger fields are observed, we believe there will arise the opportunity to accurately quantify the 
spatial structure and strength of small impact-parameter lensing events. We would expect the study of 
these events to provide a cosmological signal of small-scale substructure, and perhaps even be able to 
measure the extent to which sub-haloes are stripped. 
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